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ABSTRACT 

It has been recently shown (Tsiklauri & Viollier, 1998a) that the matter concentra- 
tion inferred from observed stellar motion at the galactic center (Eckart & Genzel, 
1997, MNRAS, 284, 576 and Genzel et al., 1996, ApJ, 472, 153) is consistent with a 
supermassive object of 2.5 x 10^ solar masses, composed of self-gravitating, degenerate 
heavy neutrinos. It has been furthermore suggested (Tsiklauri & Viollier, 1998a) that 
the neutrino ball scenario may have an advantage that it could possibly explain the 
so-called "blackness problem" of the galactic center. Here, we present a quantitative 
investigation of this statement, by calculating the emitted spectrum of Sgr A* in the 
framework of standard accretion disk theory. 
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1 INTRODUCTION 

The enigmatic radio source Sgr A* at the galactic center has 
been a longstanding puzzle. Observations of stellar motions 
at the galactic center (Eckart & Genzel, 1997; Genzel et al., 
1996) and low proper motion (< 20 km sec~^; Backer, 1996) 
of Sgr A* indicate that, on the one hand, it is a massive (2.5± 
0.4) X 10^ M0 object dominating the gravitational potential 
in the inner < 0.5 pc region of the galaxy. On the other 
hand, observations of stellar winds and other gas flows in 
the vicinity of Sgr A* suggest that the mass accretion rate 
M is about 6 x 10~^A/oyr"^ (Genzel et al., 1994). This 
implies that the luminosity of the central object should be 
more than 10*" erg sec~^, provided the radiative efficiency is 
the customary 10%. However, observations indicate that the 
bolometric luminosity is actually less than 10^^ erg sec~^. 
This discrepancy has been a source of exhaustive debate in 
the recent past. The broad-band emission spectrum of Sgr 
A* can be reproduced either in the quasi-spherical accretion 
model (Melia, 1992, 1994) with M ~ 2 x W'^Mq yr'^ 
or by a combination of disk plus radio-jet model (Falcke 
et al., 1993a, 1993b). As pointed out by Falcke and Melia 
(1997), quasi-spherical accretion seems unavoidable at large 
radii, but the low actual luminosity of Sgr A* points toward 
a much lower accretion rate in a starving disk. Therefore, 
Sgr A* can be described by a model of a fossil disk fed by 
quasi-spherical accretion. Another successful model which 
is consistent with the observed emission spectrum of Sgr 
A* has been developed by Narayan et al., 1995, 1998 and 
independently by Manmoto et al., 1997. This model is based 
on the concept of advection dominated accretion flow, in 



which most of the energy released by viscosity in the disk is 
carried along with the gas and lost into the black hole, while 
only a small fraction is actually radiated off. 

Recently, Tsiklauri & Viollier (1998a) have proposed an 
alternative model for the mass distribution at the galactic 
center in which the customary supermassive black hole is 
replaced by a ball composed of self-gravitating, degenerate 
neutrinos. It was shown that a neutrino ball with a mass 
2.5 X 10^ Mq, composed of neutrinos and antineutrinos with 
masses > 12.0 keV/c^ for g = 2 or > 14.3 keV/c^ 
for g = 1, where g is the spin degeneracy factor, is con- 
sistent with the current observational data. The purpose of 
this paper is to present calculations of the spectrum emit- 
ted by Sgr A* in the framework of standard accretion disk 
theory, assuming that Sgr A* is a neutrino ball with the 
abovementioned physical properties, and to show that this 
could resolve the " blackness problem" . 

In the recent past, Viollier et al. have proposed that 
massive, self-gravitating, degenerate neutrinos arranged in 
balls, where the degeneracy pressure balances self-gravity, 
can form long-lived conflgurations that could mimic the 
properties of dark matter at the centers of galaxies (Viollier, 
1994; Viollier et al., 1993; Viollier et al., 1992). Tsiklauri & 
Viollier (1996) have shown that a neutrino ball could play a 
similar role as a stellar cluster in the 3C 273 quasar, reveal- 
ing its presence through the infrared bump in the emitted 
spectrum. Tsiklauri & Viollier (1998b) further investigated 
the formation and time evolution of neutrino balls via two 
competing processes: annihilation of the particle-antiparticle 
pairs via weak interaction and spherical (Bondi) accretion of 
these particles. Bilic & Viollier (1997) showed how the neu- 
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trino balls could form via a first-ordor phase transition of a 
system of self-gravitating neutrinos in the presence of a large 
radiation density background, based on the Thomas-Fermi 
model at finite temperature. They find that, by cooling a 
non-dcgcncrate gas of massive neutrinos below a certain crit- 
ical temperature, a condensed phase emerges, consisting of 
quasi-degenerate supermassive neutrino balls. General relar- 
tivistic effects in the study of the gravitational phase tran- 
sition in the framework of the Thomas-Fermi model at fi- 
nite temperature were taken into account in Bilic & Viollier 
(1998a). A theorem was proven by BiUc & ViolUer (1998b) 
which in brief states that the extremization of the free en- 
ergy functional of a system of sclf-gravitating fcrmions, de- 
scribed by the general relativistic Thomas-Fermi model, is 
equivalent to solving Einstein's field equations. 



2 THE MODEL 

The basic equations which govern the structure of cold neu- 
trino balls have been derived in the series of papers (Viollier, 
1994; Violher et al., 1993; Viollier et al., 1992 and Tsiklauri 
& Viollier, 1996); here we adopt the notation of Tsiklauri 
and Viollier (1996). In this notation the enclosed mass of 
the neutrinos and antineutrinos within a radius r = of 
a neutrino ball is given by 



and fi'(r) assume the values 



M< = 87rpcr-„ 



Snpcrl {-fe') , (1) 



where 9{^) is the standard solution of the Lane-Emden equa- 
tion with polytropic index 3/2, r„ is the Lanc-Emden unit 
of length and pc is the central density of the neutrino ball. 

In the standard theory of steady and geometrically thin 
accretion disks, the power liberated in the disk per unit area 
is given by (Perry & Williams, 1993) 



D{r) 



mm 



(2) 



Here O is the angular velocity of the accreting matter, Ri is 
the inner edge of the disk and ftj denotes the angular velocity 
at the radius where its derivative with respect to r vanishes 
due to the deviation from the Keplerian law of rotation. 
Finally, the prime denotes the derivative with respect to r. 
Since the motion of accreting matter in the bulk of the disk 
is Keplerian, we assume that the angular velocity is given 
by 



fi(r) 



GM<(r) 



(3) 



In the case of a back hole M<(r) = const = Mbh, whereas 
in our case M<(r) is determined by Eq.(l). Throughout 
this paper we take the outer radius of the disk as 10^ 
Schwarzschild radii, since for larger radii the disk is unstable 
against self-gravity (e.g. Narayan et al., 1998). The radius 
of a neutrino ball with a mass 2.5 x 10® M©, composed of 
neutrinos and antineutrinos with masses niy = 12.0 keV/c^ 
for g = 2 or rriv = 14.3 keV/c^ for g — 1, is equal to 1.06 
X 10 ' Schwarzschild radii of a black hole with the same mass, 
thus the accretion disk is fully immersed in the neutrino ball. 
Moreover, as in our case there is no last stable orbit, accre- 
tion may in principle continue as r tends to zero, where n(r) 



n(0) = 



SttGpc 



a'(o) = 



(4) 



3 ' rnfi(O) ■ 

Of course the latter result is of rather academic interest, 
because in reality, the accreting matter will be diverted at 
the origin in the form of an outflow which will inevitably 
stream away perpendicular to the disk plane. The excess 
matter which has spiraled down to the very center will be 
pushed out of the plane due to the gas pressure of the ac- 
creting matter in the disk. It is important to note that this 
outflow will differ considerably from a jet shooting out of 
an accretion disk around a black hole. In the latter case, the 
jets manifest themselves as strong emitters mostly in the ra- 
dio band due to the synchrotron radiation produced by the 
electrons moving at highly relativistic velocities, whereas the 
outflow from the accretion disk immersed in a neutrino ball 
will be practically unobservable, since the outflowing mat- 
ter will be cold as it radiated off its energy while spiraling 
down in the disk (see further Fig. 2). Moreover, the parti- 
cles will be moving at non-relativistic velocities because of 
the shallowness of the gravitational potential of the neutrino 
ball that is much more spatially extended than a black hole. 
Also, it is worthwhile to note that, even at a constant accre- 
tion rate of 6 x 1O~^M0 yr~^, the baryonic mass acquired 
by the neutrino ball within the age of the universe of 10 Gyr 
would be of the order of 6 x 10* which is small compared 
to the mass of the neutrino ball. 

Numerical analysis shows that initially, as the matter 
spirals towards the center, f2'(r) is negative. From Eq.(4) 
we gather that the central value for f2' is finite and positive, 
thus there exits a point at which Q' crosses zero. This is 
precisely the point where the angular velocity attains its 
maximal value. Numerically, this happens at = i?f/r„ = 
8.25 X 10~*. Note that this position is quite close to the 
center of the ball since its radius in dimensionless units is 
Ci = 3.65375 (Cox & GiuU, 1968). Such a behavior of n(r) is 
quite interesting since, in the neutrino ball scenario, there is 
neither last stable orbit (as in the case of a black hole) nor a 
stiff stellar surface (as in the case of accretion onto a neutron 
star) . Basically, this is a consequence of the non-trivial mass 
distribution determined by the Lane-Emden equation. 

We now assume that the gravitational binding energy 
released is immediately radiated away locally according to 
the Stefan-Boltzmann's law 

D{r) = aJt^ir), (5) 

with a denoting Stcfan-Boltzmann constant. The effective 
temperature can be derived using Eqs.(l-3) and (5) yielding 

1/4 



Teff(5) 



MQQ'r„ (^g'"^ + 39') 



1 



(6) 



Here we have introduced the quantities 



n = 



2.5 X 10« MqG 



2.5 X IO'^MqC 



(-£,'^0')i = 2.71406 (Cox & Giuh, 1968) and Hi = 

Once the temperature distribution in the accretion disk 
is specified, we may calculate its luminosity using 

c2 4 exp[/»./febT(01-l' ^ ' 
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where h is Planck's constant, fcb denotes Boltzmann's con- 
stant and i is the disk inclination angle which we assume to 
be 60° as in Narayan et al. (1998). Following the same paper, 
we parameterize the accretion rate in terms of the Edding- 
ton limit accretion rate, i.e. M — rhMEddMoyr^^, where 
MEdd = lOiEdd/c^ = 1.39 X 1O^*(M/M0)g sec-^=2.21 x 
lO"^(M/M0)MGyr"\ Melia (1992) has estimated M as 
w 2 X 10~''M©yr~^ using 600 km sec~^ for the wind velocity, 
whereas Genzel et al. (1994) obtained M « 6 x lO'^Moyr"^ 
using 1000 km sec^^ for the wind velocity. These values 
translate into lO"** < m < 4 x 10~^ in terms of the Edding- 
ton units. Following again Narayan et al., 1998 we use these 
two values as the lower and upper limits for this quantity. 



3 DISCUSSION 

Results of our numerical calculations arc presented in Fig.l, 
where we plot the quantity uL^, calculated using Eq.(7). 
Data points are taken from Table 1 in Narayan et al., 1998. 
The thick solid line corresponds to the case of a neutrino ball 
with rh = 4 X 10~^, whereas the thin solid line corresponds 
to rh = 10^*. The short-dashed line represents the calcula- 
tion with a 2.5 x 10® black hole with rh = 10~^ and an 
accretion disk extending from from 3 to 10® Schwarzschild 
radii. The long-dashed line corresponds to the case when rh 
is artificially brought down to As we see from Fig.l, 

and as also was pointed out by Narayan et al., 1998, the lat- 
ter two curves provide a poor fit to the observational data. 
Actually, this is the major reason why the standard accre- 
tion disk theory was abandoned as a possible candidate for 
the description of the emitted spectrum from Sgr A*. How- 
ever, as originally was pointed out in Tsiklauri & VioUicr 
(1998a) in the neutrino ball scenario, the accreting matter 
experiences a much shallower gravitational potential than in 
the case of the black hole with the same mass, and therefore 
less viscous torque will be exerted. The radius of a neu- 
trino ball of total mass 2.5 x IQ^Mq, which is composed of 
self-gravitating, degenerate neutrinos and antineutrinos of 
mass rrii, = 12.0 keV/c^ for g' = 2 or mv = 14.3 keV/c^ for 
g — 1, is 1.06 x 10® larger than the Schwarzschild radius of 
a black hole of the same mass. In this context it is impor- 
tant to note that the accretion radius Ra = 2GM/v1, for 
the neutrino ball, where ~ 700 km/sec is the velocity 
of the wind from the IRS 16 stars, is approximately 0.02pc 
(Coker & Melia, 1997), which is slightly less than the radius 
of the neutrino ball, i.e. 0.02545 pc (for = 12.0 keV/c^ 
for g = 2 or rrii, = 14.3 keV/c^ for g = 1). Therefore, in the 
neutrino ball scenario, the captured accreting matter will 
always experience a gravitational pull from a mass less than 
the total mass of the ball. One can see from Fig. 1, that 
for this very reason the theoretical spectrum in the case of 
the neutrino ball with rh = 4 x 10~^ gives a much better fit 
than in the case of a black hole for any (even unrcalistically 
lowered) values of rh. Discrepancies between the theoretical 
and observed spectra appear in the case of the neutrino ball 
for frequencies < 40 GHz and > lO" Hz. 

At the higher end (> lO^'' Hz) of the spectrum, the 
discrepancy is due to the fact that our model does not in- 
corporate effects of Compton-scattered synchrotron radia- 
tion (which causes the second peak; on the left in Fig.l of 
Narayan et al., 1998). Our model is based on the simple- 



minded assumption of a steady, geometrically thin accretion 
disk which radiates off the gravitational binding energy lo- 
cally, according to the black-body radiation law. However, 
even in this simplified framework, our model gives a reason- 
able fit in the radio to near infrared part of the spectrum. 
Besides, it is important to note that, as it has been shown 
by Falcke & MeUa (1997), the evolution of am accretion disk 
can be considerably influenced by the deposition of mass 
and angular momentum by an infalling Bondi-Hoyle wind. 
The major result of their paper is that the modification of 
the standard accretion disk model, by taking into account 
the contribution from the Bondi-Hoyle wind and consider- 
ing the physical picture of accretion process in dynamics, 
yields significant changes in the emitted spectrum. In fact, 
it produces an infrared bump, in addition to the Big Blue 
Bump, due to the deposition of energy in the outer part of 
the fossil accretion disk. Our paper is based on the stan- 
dard accretion disk model i.e. without modifications arising 
from taking into account effects from the wind. In our case 
the gravitational potential is shallower than in the case of 
a supermassive black hole with the same mass. Therefore, 
taking into account effects from the Bondi-Hoyle wind and 
considering the non-steady problem (as in the case of Fal- 
cke & Melia's paper), both bumps will be shifted into the 
lower frequency domain. Thus the incorporation of Falcke & 
Melia's model of the accreting flow into our scenario of the 
dark matter distribution at the galactic center would pre- 
sumably produce a better fit in the < 40 GHz part of the 
spectrum. 

It is important to address the issue of consistency of 
our model with intrinsic source size versus frequency data. 
For the test we take the data of emission wavelength A = 
7 mm (Bower & Backer, 1998) and 3.5 mm (Rogers et al., 
1994; Krichbaum et al., 1994). The upper limits on the in- 
trinsic source size are < 4.1 AU (Bower & Backer, 1998) 
for 7 mm and < 1.1 AU (Rogers et al., 1994) and 2.8 ± 1.2 
AU (Krichbaum et al., 1994) for 3.5 mm assuming a dis- 
tance to the galactic center of 8.5 kpc. Now, we have to 
estimate the radial location of the circles of the accretion 
disk in our model, emitting at these two wavelengths. For 
this purpose we assume that the corresponding temperature 
of a circle can be determined by the Wien displacement law 
Vm w ?>kbT/h w 6 X 10^°T Hz (Lang, 1974), i.e. we assume 
that the maximal frequency (wavelength) in the brightness 
distribution given the black body law determines the tem- 
perature of the emitting region. This assumption seems rea- 
sonable recalling the sharpness of the maxima in the bright- 
ness distribution (see Fig. 1 in Lang, 1974). Therefore we 
obtain T7mm = Q.71K and Ts.smm = 1.437^. To find out to 
which values of the radial distance in the accretion disk these 
two values correspond, we have to use Eq.(6), graphically de- 
picted in the Fig.2. The values axe ^7mm = 8.71 x lO""* and 
^3.5mm = 1.145 X 10~*. The final predictions of our model 
would be twice these values [diameter of the emitting cir- 
cle) which in dimensional units are 2.50 AU and 3.29 AU for 
7 mm and 3.5 mm, respectively. Thus we conclude that for 
7mm our model is consistent with the observations by Bower 
& Backer (1998); the same applies to the data by Krichbaum 
et al. (1994) for 3.5 mm. However, currently some of the 
VLBI observations at millimeter wavelength stand in con- 
flict with each other (Bower & Backer, 1998). Therefore, the 
discrepancy of our estimates with the Rogers et al., (1994) 
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data is a matter of debate. Another important requirement 
which our model does satisfy is the lower limit on the size de- 
rived from the scintillation experiments (Gwinn et al., 1991). 
These experiments imply that the source diameter should be 
> 0.1 AU for 0.8 mm wavelength. Our estimates show that 
at this wavelength the source diameter is 34.11 AU which 
corroborates the validity of our model. 

As we can conclude from the latter paragraph, our 
model satisfies the source versus frequency constraints. How- 
ever, it seems unlikely that such low temperatures as re- 
quired by our model, especially close to the center of the 
neutrino ball, are actually realized, as the Galactic cen- 
ter is immersed in a hot radiation field. So far, the low- 
est termporaturc of one of the (black body radiation) emit- 
ting components introduced by Zylka, Mezger & Lesh 1992, 
was estimated to be 30 K, as a possible alternative to the 
self-absorbed synchrotron emission. Thus, we do not claim 
that in the mm wave-length band (which corresponds to 
disk temperatures of order few Kelvin) our model is capable 
of explaming the emitted spectrum reliably. The fit of the 
spectrum at these frequences (< 40 GHz), based on stan- 
dard accretion disk theory for a disk immersed in the field 
of the neutrino ball, is anyway not good, and it seems that 
some other, presumably non-thermal mechanism is respon- 
sible for the radio emission of Sgr A* . Moreover, as has been 
shown by Reynolds & McKee 1980, the compact radio source 
at the Galactic center could be a pulsar, with total emitting 
power comparable to that of the Crab pulsar. These authors 
showed that Sgr A* can be understood in the framework of a 
number of dynamically self-consistent models of incoherent 
synchrotron sources which arc energetically comparable with 
to the energy output of a few x 10'^** erg sec^^ like the Crab 
pulsar. The radio pulsar, together with a Shakura-Sunyaev 
disk embedded in the shallow gravitational potential of the 
neutrino ball, may be the key to the understanding of the 
emission spectrum of Sgr A*. In this context, it is perhaps 
important to point out that the low proper motion of the 
central radio source (< 20 km sec~^; Backer, 1996) could be 
explained by a slowly moving pulsar near the minimum of 
the gravitational potential of the neutrino ball, a possibility 
which of course would be excluded in a supermassive black 
hole scenario. 

Apart from emission of high-energy radiation by the 
pulsar, and X-ray emission of the neutrino ball through rar 
diative decay of the constituent neutrinos into light neutri- 
nos, our model would be incapable of describing of X-rays 
and gamma-rays. First, even standard accretion disk theory 
around a central black hole cannot account for the emission 
of radiation above ultra-violet. The simplest model which 
might produce X-rays is the two-component plasma model 
by Shapiro, Lightman & Eardley 1976. However, in the case 
of a neutrino ball with the physical parameters mentioned 
above, it is impossible to get X-rays and gamma-rays from 
the accretion disk by definition. In order to get an apprecia- 
ble fraction of the rest mass of the electron converted into 
X-rays, the particles need to reach a sizable fraction of the 
velocity of light, which would be impossible in our scenario 
as the escape velocity from the center of the neutrino ball 
is about 1400 km s~^. This is, in fact, the reason why our 
model can explain the blackness problem of the Galactic 
center. 

Let us now take an unprejudiced view on the high- 



energy data: In the 0.8-2.5 keV band the data available 
from ROSAT (Predehl & Triimper 1994) have a resolution 
of ~ 20 ; in the 2-10 keV band data available from ASCA 
(Koyama et al. 1996) the resolution is ~ 1'. The 35-150 kcV 
data from SIGMA (Goldwrum et al. 1994) have a resolution 
of - 15', while the EGRET data (Merck et al. 1996) from 30 
MeV to 10 GeV have a resolution of ~ 1° . The intrinsic size 
of the accretion disk (which is about the size of the neutrino 
ball) is about 0.65" which is much smaller than the reso- 
lution of current X-ray or gamma-ray detectors. Even the 
measurement by Predehl & Triimper 1994, who established 
that the X-ray source is within a 10" (0.5 pc) distance from 
the Sgr A* is not conclusive. As discussed in detail in Tsik- 
lauri & Viollier 1998a, a neutrino ball would also produce 
X-rays via the radiative decay of the heavy into light neu- 
trinos. Tsiklauri & Viollier 1998a have estimated that the 
luminosity of this emission line at energy ~ m^c^ /2 (which 
would not be a sharp lino duo to scattering by the existing 
matter at the Galactic center) should be Lj < 1.45 x 10^* 
erg/sec. In fact, this luminosity is consistent with the obser- 
vations by Predehl & Triimper 1994. 

Finally, we would like to emphasize that the idea that 
Sgr A* may be an extended object rather than a super- 
massive black hole is not new (see e.g. Haller et al., 1996; 
Sanders, 1992). To our knowledge all previous such models 
assume that the extended object is of a baryonic nature, e.g. 
a very compact stellar cluster. However, it is commonly ac- 
cepted that these models face problems with stability, and 
it has been questioned whether such clusters are long-lived 
enough, based on evaporation and collision time-scales sta- 
bility criteria (for a different point of view sec Moffat, 1997). 
Our model of Sgr A* is surprisingly simple while it satisfies 
all current observational constraints: First, a neutrino ball is 
a stable object quite alike an ordinary baryonic star, though 
much more massive, with the difference that its self-gravity 
is compensated by the degeneracy pressure of the neutri- 
nos rather than thermal pressure as in the case of a bary- 
onic star. Second, a neutrino ball with the abovementioned 
physical parameters is compact enough as to be virtually 
indistinguishable from a 2.5 x 10^ M© black hole with cur- 
rent observational resolution (~ 10^ Schwarzschild radii) of 
the observations of proper stellar motions (Eckart & Genzel, 
1997; Genzel et al., 1996). Third, a neutrino ball of this mass 
can explain its low proper motion (< 20 km sec~^; Backer, 
1996). Fourth, as a bonus of our model, the neutrino ball is 
extended enough to provide a much shallower gravitational 
potential than a 2.5 x IO^Mq black hole for the accreting 
matter, thus producing a reasonable emission flux. 
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Figure captions 

Fig.l. Comparison of tlieoretical and observed spectra 
of Sgr A*. Tlie tliick solid line corresponds to the case of a 
neutrino ball of total mass 2.5 X 10'' M© with m = 4 X 10-^ 
while the thin solid line represents rh = 10~*. The short- 
dashed line describes the calculation with a 2.5 x W^Mq 
black hole, with rh = 10~* and an accretion disk extending 
from 3 to 10^ Schwarzschild radii. The long-dashed line cor- 
responds to the case when m is artificially reduced to 10~^. 
Data points in the < 40 GHz region are upper bounds. Note, 
that the thick solid line fits the most reliable data points 
with the error bars. 

Fig. 2. Temperature of the accretion disk as a function 
of radial distance from the center. The thick line corresponds 
to the case of a neutrino ball with 2.5 x W^Mq and m = 
4 x 10"'', whereas, the thin line corresponds to a black hole 
with the same mass and accretion rate. The unit of length 
is r„ — 2.14934 x 10^^ cm. Note that for values larger than 
3.65375, which corresponds to the radius of the neutrino ball, 
the back hole and neutrino ball lines do overlap because the 
potentials are equal. This paper has been produced using the 

Royal Astronomical Society /Blackwell Science I^TJtX style 
file. 
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